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Very High Energy Cosmic Rays and Their Interactions 

Ralph Engel a 

a Forschungszentrum Karlsruhe, Institut fur Kernphysik, 76021 Karlsruhe, Germany 

The investigation of high energy cosmic rays and their interactions is a very active field of research. This article 
summarizes the progress made during the last years as reflected by the contributions to the XIII International 
Symposium on Very High Energy Cosmic Ray Interactions held in Pylos, Greece. 
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1. Introduction 

The investigation of very high energy cosmic 
rays and their interactions are inherently con- 
nected subjects of astroparticle physics. On one 
hand the understanding of cosmic ray interactions 
is needed to study the flux, acceleration and prop- 
agation of cosmic rays. For example, the high en- 
ergy cosmic ray flux can only be measured by link- 
ing secondary particle cascades observed in detec- 
tors or the Earth's atmosphere to primary parti- 
cles of certain energy, mass number and arrival 
direction. Furthermore the knowledge of parti- 
cle production is needed for the interpretation of 
secondary particle fluxes due to cosmic ray in- 
teractions in various astrophysical environments. 
On the other hand cosmic rays provide us with 
a continuous beam of high energy particles that 
can be exploited for studies of interaction physics 
at energies and phase space regions not accessible 
at man-made accelerators. 

Cosmic ray research of the last years is charac- 
terized by substantial progress in measuring pri- 
mary and secondary particle fluxes. 

Examples for new results on the primary cos- 
mic ray flux are the measurements below the knee 
by AMS, BESS and ATIC, in the knee energy 
region by KASCADE and TIBET, and at the 
highest energies by the High Resolution Fly's Eye 
(HiRes) experiment. Still some of the experimen- 
tal results appear contradictory and are subject 
of controversial discussions. For example, the re- 
sults of the composition analyses of the KAS- 
CADE and EAS-TOP data seem to be in variance 
with a first, preliminary analysis of the TIBET 
data. Similarly, there appears to be a discrep- 



ancy between the AGASA measurements of the 
cosmic ray flux above 10 20 eV and the new HiRes 
data. 

In addition to the measurement of the pri- 
mary cosmic ray flux the most powerful method 
of improving our understanding of cosmic ray 
physics is the study of secondary particle fluxes. 
New instruments measuring gamma-rays (CAN- 
GAROO, HESS, MAGIC, VERITAS, and Mila- 
gro), muons and neutrinos (AMANDA, BAIKAL, 
NESTOR, and ANTARES) have begun taking 
data or successfully performed prototype runs. 
Regarding cosmic ray physics, they are expected 
not only to test models of cosmic ray acceleration 
and interaction in supernova remnants and other 
astrophysical objects but also to provide valuable 
clues on cosmic ray composition and the charac- 
teristics of high energy particle production. 

There are many efforts to develop better mod- 
els for cosmic ray interactions or to derive in- 
formation on hadronic multiparticle production. 
The progress in this field is closely linked to mea- 
surements of forward multiparticle production in 
fixed-target and collider experiments. 

One of the central problems is the consistent 
implementation of the consequences of the steeply 
rising parton densities measured in deep inelas- 
tic e-p collisions at HERA and the indications of 
parton density saturation seen at the Relativis- 
tic Heavy Ion Collider RHIC. RHIC data clearly 
demonstrate the difficulties of extrapolating mod- 
els tuned to accelerator data to higher energy or 
other projectile/target combinations. Many mod- 
els predicted a secondary particle multiplicity ex- 
ceeding that measured in central Au-Au collisions 
by ~ 30% or more. The impact of the RHIC 
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data on the extrapolation of cosmic ray interac- 
tion models to ultra-high energy is still far from 
being understood. 

Measurements of cosmic ray showers and sec- 
ondary particle fluxes have reached a preci- 
sion that they become increasingly important in 
constraining hadronic interaction models. De- 
spite providing mainly indirect information on 
hadronic multiparticle production they allow the 
exclusion of extreme model extrapolations and 
limit exotic physics scenarios. 

This article presents a summary of recent de- 
velopments in the field of very high energy cos- 
mic ray physics and related interaction physics, 
focusing on the contributions presented at the 
XIII International Symposium on Very High En- 
ergy Cosmic Ray Interactions. The plan of the 
paper is as follows. In Sec. El the current status 
of cosmic ray flux measurements is given. Re- 
sults of different measurements are compared and 
their dependence on hadronic interaction mod- 
els employed for data analysis is discussed. The 
progress in modeling extensive air showers is out- 
lined in Sec. focussing on status and uncer- 
tainties of high-energy interaction models. Moti- 
vated by the current use of QGS JET as "standard 
candle" interaction model in almost all high en- 
ergy cosmic ray experiments, uncertainties and 
features of interaction models are discussed in 
some detail. The importance of analyzing ob- 
servables of relevance to cosmic ray physics in 
experiments at current and future accelerators is 
emphasized. Sec. 0] gives a short update on the 
controversial subject of exotic interaction features 
claimed to be found in emulsion chamber mea- 
surements. The very active field of measuring 
gamma-rays, muons and neutrinos produced in 
cosmic ray interactions is briefly touched upon 
in Sec. El Finally, conclusions and an outlook is 
given in Sec. El 

2. Cosmic ray flux 

For understanding very high energy cosmic rays 
and their sources, the measurement and inter- 
pretation of the all-particle flux, the elemental 
composition (including the 7-ray fraction) , the ar- 
rival direction distribution (large scale anisotropy, 



small scale clustering, and correlation with hypo- 
thetical sources), and temporal variations are of 
central importance. We will briefly discuss these 
topics beginning with the highest energies. Re- 
cent reviews of the experimental situation can be 
found in HEEEEI. 

2.1. Ultra-high energy cosmic rays 

At the highest energies, the Akeno Giant Air 
Shower Array (AGASA) and the High Resolution 
Fly's Eye (HiRes) detectors are the installations 
with the biggest accumulated aperture that have 
published flux data. The flux data from both 
experiments are shown in Fig. ^ 

First of all there is the well-known and often 
discussed discrepancy between the two data sets 
at very high energy. Whereas the AGASA data 
do not show any sign of the expected GZK cut- 
off [21] in the energy spectrum |7I25| . the HiRes 
results are compatible with such a GZK signa- 
ture |10lll| . The statistical significance of this 
discrepancy above 10 20 eV is about 2-3 a At 
lower energy, the overall difference between the 
measurements is well within the range of the sys- 
tematic errors of both experiments |27l28j . This 
also applies to the independent data set from 
the Yakutsk array |2H] which is characterized by 
a larger shower-to-shower reconstruction uncer- 
tainty of 32 - 46% as compared to about 20% 
for HiRes and AGASA. The Yakutsk array also 
has an integrated aperture ~ 40% smaller than 
AGASA. 

Secondly there is a smooth transition between 
the HiRes and AGASA data and the lower en- 
ergy measurements of the prototype instrument 
HiRes-MIA [EJ and the Akeno air shower array 
|8I9| . respectively. This might indicate a system- 
atic bias in one or both measurement techniques 
that could be related to the simulation of ultra- 
high energy air showers. 

To estimate the elemental composition of ultra- 
high energy cosmic rays (UHECR), both AGASA 
and HiRes have analyzed their data in terms of a 
two-component proton/iron composition hypoth- 
esis. 

The HiRes analysis is based on the measure- 
ment of the depth of shower maximum. Using 
the hadronic interaction model QGSJET |30I31| 
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Figure 1. Primary cosmic ray flux scaled with E 2,5 . Shown is a selection of recent measurements 
as discussed at this meeting together with some older data for reference (AGASA 0, Akeno |7I8I9| . 
HiRes HiRes-MIA [B], KASCADE fHIT5| . MSU QU, RUNJOB fT71T3| . ATIC OH). For the 

sake of clarity, the all-particle fluxes from EAS-TOP [20121) and Tibet [22123) are not shown. They 
cannot be distinguished from the others in this representation. 



for interpreting the data they find a transition 
to a light composition |13| that remains proton 
dominated (80% protons) between 10 17 eV and 
10 i9.3 e y [K)!;^) . Preliminary results of the re- 
analysis of the AGASA muon density data mea- 
sured with the Akeno muon detectors show also 
a transition to a proton dominated composition. 
In contrast to the HiRes results, the transition 
is found to occur gradually over a large energy 
range. From about 50% iron fraction at 10 175 eV 
the iron contribution drops below 20% at 10 19 eV 
|7I33| . Similar to the HiRes analysis, the Yakutsk 
composition measurement is related to the 
primary mass sensitivity of the shower depth of 
maximum. The Yakutsk group find a light com- 
position of about 70-80% proton and helium in 
the energy range 5 x 10 17 - 5 x 10 18 eV 

The most natural interpretation of the chang- 



ing composition would be the transition from 
Galactic to extragalactic cosmic rays. The tran- 
sition energy would be somewhere between 10 17 
and 10 19 eV. Whereas HiRes data are consistent 
with the interpretation that the "ankle" in the 
cosmic ray spectrum is already a signature of 
the GZK cutoff (e + e~-pair production of protons 
with photons of the CMB) [35136) . AGASA data 
favour an interpretation of the ankle as transition 
region between Galactic and extragalactic cosmic 
rays. 

It is clear, however, that these composition 
measurements are strongly model dependent as 
there is a large theoretical uncertainty in pre- 
dicting electron and muon shower sizes as well as 
the depth of shower maximum for hadron- induced 
showers 37 38 39 . An interpretation based on 
SIBYLL [40141) or neXus gives a heavier com- 
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position: about 30 and 50% iron, respectively (see 
also Fig. [5J. It is unclear which of the model 
predictions is more realistic and also whether the 
range of predictions exhausts the range of the the- 
oretical uncertainties (see discussion in Sec. |2J. 
Moreover there are signs of inconsistencies |43I3| . 
For example, the muon densities measured for the 
same showers as used in the depth of maximum 
analysis of the HiRes-MIA data are similar to 
or even exceed those expected for iron primaries 
[T5]. Furthermore investigations based on mass- 
sensitive observables such as shower disk thick- 
ness and shape of the lateral distribution also in- 
dicate a heavier composition (~ 80 — 90% iron) 

Limits on the primary photon fraction were 
given by AGASA based on the investigation of 
showers with muon density information. The 
fraction of photon-induced showers is found to be 
smaller than 28% (67%) at energies greater than 
10 i9 ( 10 i9-5) eV at 95% CL. A recently developed 
method of comparing shower-by-shower measure- 
ments with theoretical expectations for photon- 
induced showers |46| allows one to derive a limit 
at even higher energy where the statistics is very 
sparse: less than 65% of all showers with ener- 
gies above 1.25 x 10 20 eV are initiated by photons 
(95% CL) HZJ. 

Given the limited statistics accumulated until 
now, the arrival direction distribution of UHECR 
appears isotropic. There are a number of cosmic 
rays forming arrival direction multiplets in the 
AGASA data set (57 events with E > 4 x 10 19 eV: 
6 doublets, 1 triplet) |7l48j . The statistical sig- 
nificance of this small scale clustering is subject 
to controversial discussion. If the clustering were 
found with "a priori" chosen values for energy 
threshold and separation angle, i.e. without per- 
forming a scan in energy threshold and separation 
angle, the chance probability would be less than 
10 -4 . Assuming that such a scan was performed, 
the chance probability would increase to about 
0.3% ng . 

The exposure of HiRes in stereo mode has not 
yet reached that of AGASA 1 . There are 27 events 

1 Viewing showers in monoscopic mode HiRes I has reached 
a higher accumulated aperture than AGASA for energies 
above 3 X 10 19 eV. However, due to the limited and highly 



detected in stereoscopic mode above 4 x 10 19 eV. 
Using this data set in an autocorrelation analy- 
sis no significant small-scale clustering is found. 
Adding the HiRes stereo data set to that from 
AGASA only one additional pair is found. The 
clustering found in the combined data set is esti- 
mated to have a chance probability of 1%. 

The search for correlations with astrophysical 
sources is hampered by the incompleteness of cat- 
alogs and related object detection and selection 
biases. Assuming that the source distribution fol- 
lows that of the luminous matter in the universe it 
is natural to expect a correlation with the super- 
galactic plane [S3j- No such correlation is found 
in the current data sets. There are a number 
of correlations claimed between the AGASA and 
Yakutsk data sets and BL Lacs |5 1152153] . With 
an angular resolution of about 0.7°, the HiRes 
stereo data set is ideally suited for such studies. 
Using the same astrophysical objects no correla- 
tion with the HiRes events is found for energies 
above 2.4 and 4 x 10 19 eV. 

Recently new indications of a correlation with 
BL Lacs were found if the energy threshold for 
comparison is lowered to 10 19 eV (10 out of 271 
showers have arrival directions within 0.8° of one 
of the 156 selected BL Lacs from the Veron cata- 
log) The chance probability of finding such 
a correlation in an un-correlated data set is es- 
timated as 0.1%, but only a new, independent 
data set will allow to assess the significance un- 
ambiguously. Taken at face value, the correlation 
would have to be interpreted as a small fraction 
of neutral particles in UHECR. 

Currently there are two large- aperture detec- 
tors for UHECR in construction, the Pierre 
Auger Observatory |55I56| and the Telescope Ar- 
ray (TA) |57I58| . Both detector concepts em- 
ploy the hybrid detection technique of measuring 
air showers with surface detectors (Auger: water 
Cherenkov tanks, TA: plastic scintillators) and 
fluorescence telescopes. The hybrid technique 
will allow a good energy calibration of showers 
measured with surface detectors and improve the 
ability of composition measurements. 



asymmetric angular resolution, the HiRes I mono data set 
is not suited for studying small scale clustering. 
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After testing the detector design with an engi- 
neering array |59) the construction of the south- 
ern Auger observatory in Malargue, Argentina is 
in full progress At the time of writing this 
article about 800 of the planned 1600 surface de- 
tector stations and 50% of the fluorescence tele- 
scopes are completed. Already during the con- 
struction phase the integrated aperture of Auger 
has reached that of 10 years of data taking with 
AGASA. It is planned to build a similar observa- 
tory in the northern hemisphere to obtain nearly 
uniform full sky coverage. 

2.2. Knee energy region 

In contrast to previous measurements in the 
knee energy region (for example, see the compila- 
tions in 60 2 (JJ) the recent all-particle flux mea- 
surements by KASCADE Q3|, EAS-TOP |U], 
and Tibet [231 agree within 10% with each other. 
The knee is at about 3 x 10 15 eV and no deviation 
from a broken power law with a smooth transition 
region is found within the current experimental 
resolution. 

Concerning the elemental composition, the sit- 
uation is less clear. There is increasing evidence 
for a transition from a mixed to a heavy com- 
position with increasing energy. However, the 
detailed change of composition through the knee 
energy range is still unclear. 

All experimental results discussed at the 
meeting show a trend towards a heavy 
composition with increasing energy (KAS- 
CADE [H], Tibet E2, EAS-TOP [H| , and 
SPASE/ AMANDA |H21 ) - However, it is diffi- 
cult to find observables that demonstrate this 
transition beyond doubt as all composition stud- 
ies depend very much on the hadronic interaction 
models used for data interpretation and alter- 
native explanations might be possible, although 
unlikely (see, for example, Probably the 

least model dependent analysis is that of muon- 
rich and muon-poor showers by the KASCADE 
Collab. |SU, which demonstrates that the knee 
is mainly caused by a disappearance of the light 
(i.e. muon-poor) flux components. 

The high-statistics data of multi-detector setup 
KASCADE [SSj allow the analysis of the corre- 
lated muon (Efj, > 240 MeV) and electron sizes 
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Figure 2. Flux derived for 5 elemental groups 
from KASCADE data [TITTS] . The top panel 
shows the results obtained with QGSJET 01, the 
bottom panel those with SIBYLL 2.1. 



of showers in terms of 5 mass groups [14115) . 
Fig. El shows the results for two hadronic interac- 
tion models, QGSJET and SIBYLL. In this anal- 
ysis, the derived all-particle flux turns out to be 
almost independent of the used hadronic inter- 
action model. However, the different elemental 
fluxes vary strongly. For example, using QGSJET 
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the flux appears dominated by helium below the 
knee with no significant iron contribution. The 
SIBYLL-based interpretation favours instead he- 
lium and carbon below the knee and a small but 
significant fraction of iron primaries is needed. 
Both models do not provide a fully consistent 
description of the KASCADE N e -Nf, data. The 
found deviations underline the high statistical ac- 
curacy of the KASCADE data and show the need 
of improving hadronic interaction models. 

The interpretation of the KASCADE data with 
both models shows that the flux of the light com- 
ponents exhibits a break in the power law at 
different energies with lighter elements having a 
lower break energy. No spectral break is found 
for iron in the considered energy range. One 
of the central questions is that of the scaling of 
the break energies. Acceleration and propaga- 
tion models for the knee typically predict rigidity- 
dependent scaling whereas models with new par- 
ticle physics lead to mass-dependent scaling. Un- 
fortunately, the strong had. interaction model de- 
pendence does not allow us to draw conclusions 
on a possible mass- or rigidity-dependent scaling 
of the break energies. 

The composition measurement by the EAS- 
TOP Collab. |20I21 | is based on 3 elemental 
groups and uses electron and GeV-muon data. 
It shows the same qualitative behaviour of the 
elemental groups as seen in the KASCADE anal- 
ysis. The iron-like mass group does not exhibit a 
break in the power spectrum and seems to have 
a harder power law index than the light compo- 
nent. The knee appears to be caused by elements 
in the mass range of proton and helium with a 
break energy of about 3.5 x 10 15 eV. 

Some of the air showers measured with the 
EAS-TOP and SPASE arrays produce high- 
energy muon bundles that can be detected with 
MACRO (En > 1.3 TeV) and AMANDA (E^ > 
300 GeV), respectively. Again, analyses of the 
coincidence data sets show a transition from a 
mixed to a heavy composition 20 66 I62I67| . 

The preliminary and statistically limited mea- 
surement of the proton flux by the Tibet AS7 
Collab. seems to be at variance with KAS- 
CADE and EAS-TOP results. Over the energy 
range from 2 x 10 14 to 10 16 eV the proton spec- 



trum is found to follow a power law with the index 
3.14±0.10. By comparing the Tibet proton spec- 
trum with that of RUNJOB and JACEE a much 
lower break energy of ~ 5 x 10 14 eV is inferred. 

The reasons for the discrepancy between Tibet 
and KASCADE/EAS-TOP data is not yet un- 
derstood. However, it should be noted that both 
KASCADE and EAS-TOP employ in their anal- 
ysis the electron-muon size correlations in show- 
ers whereas the Tibet measurement is based on a 
neural net analysis of a number of shower observ- 
ables that combine emulsion chamber information 
with scintillator data: N 7 (multiplicity of a fam- 
ily), ^Ej (energy sum of a family), (i? 7 ) (mean 
lateral spread of a family), (Ej ■ i? 7 ) (mean lat- 
eral spread of the family energy) , and N e (shower 
size) |68| . Furthermore all three experiments are 
at different altitudes and probe different stages 
of shower evolution. Tibet is at a vertical depth 
of 606 g/cm 2 , EAS-TOP at 820 g/cm 2 and KAS- 
CADE at 1020 g/cm 2 . The use of different air 
shower simulations can be ruled out as all exper- 
iments now apply CORSIKA 

Due to the large statistical errors, the direct 
flux measurements of individual mass groups by 
JACEE 70 1, RUNJOB [T7IT5] and ATIC |7T] 
do not yet impose strong constraints on the air 
shower data. All composition data discussed here 
are compatible with possible extrapolations of di- 
rect measurements at lower energy. 

In cosmic ray models that explain the knee by 
propagation effects (leakage from our Galaxy), an 
increasing dipole anisotropy of the shower arrival 
directions is expected (i.e. [12]). Analyzing 2x 10 7 
showers in the knee energy range the KASCADE 
group do not find any significant anisotropy signal 
|73j. Similarly, no cosmic ray point sources are 
seen |73) . 

Being located at higher altitude, the Tibet ar- 
ray has a much lower energy threshold. Again, 
no Galactic anisotropy was found in the Tibet 
AS7 data. However, using more than 5 x 10 9 
showers with E > 3 x 10 12 eV, the Tibet AS7 
Collab. could detect the dipole anisotropy due to 
the orbital motion of the Earth around the Sun 
(Compton-Getting effect) with an amplitude of 
about 0.1% |22T75] (see also |7J]). 
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3. Modeling of cosmic ray interactions and 
EAS 

Given the dependence of the cosmic ray flux 
and composition measurements on the under- 
standing and modeling of hadronic interactions 
of cosmic rays and their secondary particles, it 
is natural to assume that discrepant results dis- 
cussed in the previous section can be, at least to 
some extent, traced back to the use of different 
models for data interpretation. 

Several groups have shown that an analysis of 
the same data set with different hadronic interac- 
tion models can lead to a wide range of different 
results (see, for example, |21I15| 1. This means 
that the use of the same shower simulation model 
is pre-requisite for a fair comparison of the results 
of different experiments. 

Another important aspect of inter-experiment 
comparison is the use of sufficiently realistic and 
accurate shower simulations. For example, using 
the same model for shower evolution, one can ob- 
tain different interpretations of the data if differ- 
ent observables of the showers are considered jZJ] . 
Therefore experiments might arrive at contradict- 
ing conclusions even if the same shower simula- 
tion tools are used. 

The largest uncertainty in EAS simulation 
stems from the unknown characteristics of 
hadronic multiparticle production |78l79j . As has 
been realized during the last years, also interac- 
tions at intermediate energies can contribute sig- 
nificantly to this uncertainty, though to a smaller 
extent {80181182] . In addition there are uncertain- 
ties coming from the treatment of electromagnetic 
interactions and differences in details of particle 
transport and decay implementations |79| . 

Motivated by different models of hadronic mul- 
tiparticle production, three energy regions are 
distinguished. At very low energy (close to the 
particle production threshold up to a few GeV) 
particle production is characterized by the pro- 
duction and decay of resonances. Knowing all 
resonances and their decay branching ratios al- 
lows one to construct a rather complete model 
for the interaction cross section and hadronic fi- 
nal states. For example, the codes HADRIN [53] 
and SOPHIA [SI] have many resonance chan- 



nels tabulated. The intermediate energy range 
up to about 10 3 GeV can be well understood 
in a model that describes particle production 
on the basis of the fragmentation of two color 
strings (i.e. older versions of FLUKA 85 ). At 
energies above 10 3 GeV minijet production and 
multiple parton-parton interactions become im- 
portant and require again a different modeling. 
The most frequently used high energy models are 
DPMJET |Sjj], neXus @2, QGSJET |3T], and 
SIBYLL gU|. 

In the following we will summarize some impor- 
tant recent developments in modeling hadronic 
interactions and related activities of measuring 
hadron production in accelerator experiments, 
and will discuss new trends in air shower simu- 
lation. 

3.1. Hadronic interactions at high energy 

There are basically four central assumptions 
that characterize a model's high energy extrap- 
olation |88I89I9()| : 

(i) size and energy dependence of the QCD 
minijet cross section, 

(ii) distribution of partons in transverse space 
(profile function), 

(iii) scaling of leading particle distributions or 
scaling violation, and 

(iv) treatment of nuclear effects (semi- 
superposition model, Gribov-Glauber ap- 
proximation, increased parton shadowing, 
etc.) 

It is beyond the scope of this article to discuss 
all models and their differences. We shall concen- 
trate here on aspects relevant to p-air interactions 
and consider only the two most frequently applied 
models, QGSJET and SIBYLL. Apart from the 
treatment of nucleus- nucleus collisions, QGSJET 
and SIBYLL differ mainly in the first two points 

Since version 2.1, "post-HERA" parton densi- 
ties are used in SIBYLL for calculating the mini- 
jet cross section whereas QGSJET 01 (and the 
earlier version QGSJET 98) was developed with 
older, "pre-HERA" parton densities. Another 
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Figure 3. Electron-muon size correlation for showers simulated with CORSIKA 87 . Predictions of the 
old and new versions of QGSJET jSH] are compared with SIBYLL 2.1 |41j . 



important difference between the models is the 
treatment of the minijet transverse momentum 
cutoff needed to restrict the minijet calculation 
to the perturbative QCD domain. In QGSJET 01 
an energy- independent, constant value of 2GeV 
is used. The SIBYLL authors implemented an 
energy-dependent transverse momentum cutoff 
whose value is similar to that of QGSJET 01 
at low energy but increases to about 8GeV at 
10 20 eV gT]. This was needed as "post-HERA" 
parton density functions predict at ultra-high en- 
ergy gluon densities that lead to overlap of in- 
dividual gluon wave functions in a proton for a 
transverse momentum cutoff as low as 1.5 GeV 
(see a recent review on this subject is [93 ). 
In this phase space region non-linear evolution 
equations have to be used to describe parton den- 
sities. The expected size of the non-linear correc- 
tions is theoretically not understood and subject 
of intense research There are models that 

predict an early and total saturation of the gluon 
density (color glass condensate [HSJ ) and others 
with moderate changes. 

Many experiments are searching for signs of 
deviations from linear parton density evolution 
equations or gluon density saturation. Although 



HERA and RHIC are colliders with CMS ener- 
gies of y/s ~ 200 GeV, corresponding to only 
2 x 10 13 eV, they are currently the best instru- 
ments for studying saturation effects. At HERA, 
parton densities of quarks are measured directly 
and, through scaling violations, that of gluons de- 
rived. HERA data can be described assuming 
parton density saturation, for example, in terms 
of the Golec-Biernat-Wusthoff model or ap- 
plying perturbative QCD without any non-linear 
effects 97 98 . At RHIC, parton densities can- 
not be measured directly. However, the scaling of 
jet rates and other observables with the number 
of participating nucleons depends on the assump- 
tions on the number of partons in the very gluon- 
dense environment of a heavy nucleus. Many 
aspects of RHIC data indicate strong deviations 
from naive parton model predictions 9 9ll00j but 
the energy range of the collisions is too limited to 
show unambiguously that the effects observed so 
far are requiring parton density saturation |l()lj . 

The problem of using "post-HERA" parton 
densities for extrapolating hadronic interactions 
to ultra- high energy y/s ~ 500 TeV within the 
Quark-Gluon Strings Model |102| is addressed in 
the new version of QGSJET, called QGSJET- 
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II |103| . Different from the SIBYLL approach, 
the transverse momentum cutoff is kept energy- 
independent. This is achieved by introducing 
non-linear effects for partons below the pertur- 
bative scale, equivalent to non-linear evolution 
equations. As one cannot speak of individual 
partons in the soft, non-perturbative domain, 
these non-linear effects are implemented as multi- 
pomeron interactions (enhanced pomeron graphs 
104 ), which are summed to all orders. Other 
important improvements are the treatment of 
diffraction dissociation. Whereas the old ver- 
sion of QGSJET had a fixed ratio of diffractive 
to elastic cross sections (grey disk limit, see also 
|105| L the new version approaches at high en- 
ergy the black disk limit (see discussion in |88|1. 
Furthermore QGSJET-II was tuned to better de- 
scribe the secondary particle multiplicity at low 
collision energy. Although the latter changes are 
more of technical nature they still are important. 
The old version predicted at low energy too high 
a pion multiplicity, a possible reason for higher 
GeV-muon multiplicities obtained for EAS than 
found with other models |10fll82j . It is clear that 
QGSJET-II is a much more theoretically consis- 
tent model than the previous versions. 

There are a number of important consequences 
from these changes in QGSJET. First of all the 
low- a; extrapolation of the parton densities be- 
comes less steep than naively expected in linear 
perturbative QCD. Secondly there are changes of 
the effective parton densities acting in hadron- 
hadron collisions in dependence of the projectile 
and target mass number A (suppression for large 
A). This leads to the violation of the superpo- 
sition approximation even for fully inclusive ob- 
servables. 2 Thirdly the fluctuations in inelastic- 
ity are considerably reduced. The leading parti- 
cle distributions are now qualitatively similar to 
those of SIBYLL and DPMJET [TTTTj . 

The impact of the QGSJET improvements on 
air shower predictions is currently under inves- 
tigation 107 . For example, Fig. |3] shows the 
electron-muon size correlation for vertical EAS 

2 In the simplest version of the superposition model, an 
iron-induced shower is equivalent to 56 proton induced 
showers having 1/56 of the shower energy. The superposi- 
tion model is expected to be valid for inclusive observables. 



at sea level |S7]. The contours are iso-lines of 
the correlation function at half maximum. For 
a given energy the number of muons is reduced 
significantly. At the same time 30-40% more elec- 
trons are expected at detector level in the knee en- 
ergy region. The interpretation of data from ex- 
periments that utilize electron and muon numbers 
as composition sensitive observable (e.g. KAS- 
CADE, EAS-TOP) will change towards a heavier 
composition if QGSJET-II is used. 

At the highest energies, the predictions of the 
electron numbers of the old and new QGSJET 
versions are not very different. Still it can be 
expected that the energy calibration of exper- 
iments like AGASA would have to be revised 
downward though detailed simulations are needed 
to estimate by what amount. Experiments like 
Auger |55| will be much more sensitive to the 
modifications: the predicted muon density at 600 
to 1000 m from the shower core is reduced by 
- 30%. Using QGSJET-II also the interpreta- 
tion of fluorescence measurements will change and 
shift towards a heavier composition as the mean 
depth of maximum is increased by 10g/cm 2 and 
20g/cm 2 at 10 20 eV for proton and iron showers, 
respectively [3"5] . 

Different assumptions on the QCD minijet 
cross section, i.e. calculated with or without satu- 
ration, lead to enormous differences in the model 
extrapolations [EE|- The most striking example 
is the secondary particle multiplicity in p-air col- 
lisions. At ultra-high energy, QGSJET 01 pre- 
dicts a more than 3 times higher secondary parti- 
cle multiplicity than SIBYLL |78I89| . Due to the 
implementation of new parton densities, the mul- 
tiplicity of QGSJET-II is even higher than that of 
QGSJET 01 at ultra-high energy. However, these 
differences are of very much reduced importance 
for the evolution of air showers as most of the 
secondary particles have a very small energy. Of 
greater importance are two model characteristics 
that are indirectly linked to the minijet cross sec- 
tion: the total p-air and 7r-air cross sections and 
the distribution of leading secondary particles. 

The total cross section of a model depends at 
high energy mainly on the minijet cross section 
and the transverse profile function which is a mea- 
sure of how partons are distributed in a hadron 
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Figure 4. Compilation of p-air production cross 
sections and model predictions (from [7^1, modi- 
fied and updated |1U7| ^ . 



|89|. It is the difference of the assumed profile 
functions and minijet cross sections and the lack 
of data to distinguish between these assumptions 
that lead to widely varying high-energy cross sec- 
tion extrapolations (see Fig. HJ. Cross section 
measurements over a wide range of energy would 
help to reduce the model ambiguities. Although 
hadron-air cross sections are needed for simula- 
tion, p-p and p-p cross section measurements are 
also of great interest. Using the Gribov-Glauber 
approximation, cross sections with air can be es- 
timated using nucleon-nucleon cross section data. 
At the moment, due to the contradicting total p-p 
cross section measurements at Tevatron 108 109 , 
possible theoretical extrapolations are experimen- 
tally not very much restricted |88lll()j . Therefore 
the planned total cross section measurement with 
the TOTEM/CMS detector combination at the 
LHC |111I112| . corresponding to an equivalent 
energy of about 10 17 eV, will be of outstanding 
importance. 

All current high energy interaction models are 
based on the implicit assumption that leading 
particle distributions scale with energy The lead- 



ing particle distributions are tuned to low-energy 
data and change at high energy only due to 
energy-momentum conservation effects as the en- 
ergy has to be shared between the leading par- 
ticles and the increasing bulk of low-energy sec- 
ondary particles. This assumption seems to de- 
scribe the very sparse data we have on lead- 
ing particle production up to HERA energy (~ 
2 x 10 13 eV) |113j . However, there are theoretical 
arguments that the leading particle distributions 
will change drastically at ultra- high energy |114| . 

If parton density saturation indeed occurs at 
cosmic ray energies, a collision can be viewed 
as black disk scattering: the gluons completely 
"fill" the target nucleus |37I115I116| . Not only 
the very numerous partons at small x but also 
the much faster valence quarks will participate in 
the interaction. Indeed, in a non-peripheral colli- 
sion (complete saturation) the chance probability 
of valence quarks to scatter off these gluons ap- 
proaches unity. As a consequence this will lead to 
the disintegration of the leading valence di-quark: 
no leading baryon is produced and the elasticity 
of the collision drops by almost a factor 2. 

Indeed there are some indications of "anoma- 
lous" baryon stopping in heavy ion collisions, 
however, at much lower energy (for example, |117l 
I118| ) and different theoretical interpretations are 
put forward (for example, string junction inter- 
pretation: |119I120| and saturation |121| ). Mea- 
surements of hadron production in the very for- 
ward direction at RHIC H|, Tevatron [TUB], and 
LHC |111| will be needed to study the leading 
baryon distributions systematically and clarify 
the situation. 

The implementation of different scenarios of 
parton density saturation in the SIBYLL 2.1 code 
allows a first estimate of the expected effects and 
their dependence on model parameters |115lllf>| . 
In a conservative scenario the mean depth of 
shower maximum, (X max ), is reduced by about 
20 g/cm 2 at 10 19 eV, corresponding to the dif- 
ference between SIBYLL and QGSJET 01 pre- 
dictions. Depending on the assumptions, much 
larger reductions of (X max ) are possible. 
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3.2. Hadronic interactions at intermediate 
energy 

Hadronic interactions in air showers at inter- 
mediate energy are often simulated with mod- 
els like GHEISHA Q22I (used in CORSIKA JgQ) 
or the Hillas splitting algorithm |123| (used in 
MOCCA [El] and AIRES Q2H1). Both models 
are very fast but rather crude parametrizations 
of low-energy data or interaction physics. Their 
application is certainly justified if only the elec- 
tron/photon component of a shower or calorimet- 
ric quantities are studied as in this case the de- 
tails of low-energy interactions are of minor im- 
portance. 

The situation is, however, different for muons. 
Due to successive hadronic interactions, the num- 
ber of pions and kaons increases in an air shower 
with decreasing particle energy. Below some 100 
(500) GeV pions (kaons) are more likely to de- 
cay than to undergo further interactions. There- 
fore, hadronic interactions in the energy range 
from several GeV to a few hundred GeV are very 
important for understanding GeV muon produc- 
tion in EAS of all energies |80I81| . Recent stud- 
ies have shown that the muon density at large 
lateral distance is indeed very sensitive to the 
model used for low- and intermediate-energy in- 
teractions. The differences between the predic- 
tions of the various models are of the order of 
10-20% in the relevant lateral distance range but 
can be as large as 50% [126182) . 

A detailed comparison of low- and 
intermediate-energy models to available data 
[8211 2 7j show that GHEISHA does not provide 
an adequate paramctrization of the interaction 
characteristics. The Hillas splitting algorithm 
seems to give a somewhat better description but 
a thorough comparison to data is hampered by 
the limitation to only p/-7r-air collisions. The 
best models available are clearly FLUKA [HSj 
and UrQMD |128| but there are still significant 
differences between the predictions of these two 
models. 

The energy range up to 400 GeV is in reach 
of fixed-target accelerator experiments. Not only 
that such experiments can easily measure with 
light, air- like nuclei as targets they can also run 
with tagged pion and kaon beams and measure 



particle production in the very forward direction. 
Recognizing the importance of low-energy mea- 
surements for atmospheric neutrino flux predic- 
tions |129) and neutrino factories, a programme 
was begun to systematically measure pion and 
kaon production in minimum bias collisions. Ex- 
amples are the HARP experiment |13()ll31j where 
first data are available now, the NA49 minimum 
bias p-C run |130| . and the MIPP experiment 
|108I132| which is currently taking data. 

3.3. Information from air shower measure- 
ments 

It is difficult to obtain information on hadronic 
multiparticle production at ultra-high energy 
from EAS measurements. First of all the primary 
particle mass is, in general, not known. Secondly, 
the large number of successive hadronic interac- 
tions smears out any striking features of the pri- 
mary interaction. Therefore most analyses of air 
shower data in terms of interaction physics are 
highly indirect and often serve only the exclusion 
of extreme scenarios (for example, |37|). 

One of the most interesting measurements of 
this kind is the analysis of deeply penetrating air 
showers to obtain the inelastic proton-air cross 
section. Traditionally an exponential function 
is fitted to the observed A max distribution and 
the derived absorption length A is converted to 
the interaction length via a so-called fc-factor 
(for a recent discussion, see |133| ). The HiRes 
Collab. developed a new method that does par- 
tially avoid the ambiguities of the definition of 
a fc-factor jl84j . Applying this method to the 
stereo data set, a preliminary p-air cross section 
of 456 ± 17(stat) + 39 - ll(sys)mb at 10 18 5 eV 
is derived |135| . This cross section is lower than 
the current model extrapolations, see Fig. A 
number of possible biases still need to be inves- 
tigated. For example, a small fraction of pho- 
ton primaries would be enough to spoil the cross 
section measurement as photon-initiated showers 
have a much larger depth of maximum. Further- 
more the self-consistency of the method should be 
checked by applying it to a model that is modified 
to match the actually measured cross section. 

Experiments that measure many observables of 
air showers simultaneously can check the con- 
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sistency of EAS simulation. For example, the 
KASCADE installation allows the measurement 
of shower size (electrons), muon densities with 
thresholds of 240, 490, and 2400 MeV, and hadron 
multiplicities and energies above 70 GeV in the 
shower core [141136) . The correlation between the 
different observables provides constraints on in- 
teraction models even if the full range of possible 
primary particles is considered [137) . The latest 
versions of the had. interaction models available 
in CORSIKA satisfy these constraints but some 
earlier versions are clearly at variance with the 
data. 

Emulsion chamber experiments at high altitude 
have the advantage that they have a low shower 
energy threshold, reaching into the region of di- 
rect primary flux measurements. Therefore they 
are well suited for testing shower simulation mod- 
els. For example, by comparing the predicted and 
measured optical density distribution (i.e. the en- 
ergy distribution of particles in the TeV range) in 
emulsions of the Pamir experiment, it was found 
that some old versions of hadronic interaction 
models could be excluded (!H8llH9j . 

Finally it should be noted that there are some 
indications of a systematic discrepancy between 
the mass composition derived from N e — iV M based 
measurements and that from data sensitive to 
the depth of maximum [1401141] . It is unclear 
to what extent these discrepancies are related to 
shortcomings in the simulation of hadronic inter- 
actions, but one can try to bring different mea- 
surements into better agreement by modifying the 
underlying shower simulation correspondingly. A 
hadronic interaction model with small cross sec- 
tion and a somewhat reduced inelasticity as com- 
pared to QGSJET 01 is favoured in this analysis 

EH- 

3.4. Trends in EAS simulation 

There is a lack of statistics of Monte Carlo sim- 
ulations in many comparisons of data with theo- 
retical predictions. For example, the statistics of 
the KASCADE shower data exceed the simulated 
one available for each interaction model combi- 
nation by about a factor 10. A similar disparity 
of data to simulation statistics is also typical for 
modern large-scale detectors such as Auger. The 
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Figure 5. Mean depth of shower maximum for 
proton- and iron-induced showers [142| . The pre- 
dictions of different hadronic interaction models 
as calculated with different shower codes are com- 
pared. 



problem is even amplified if many observables are 
used to characterize each shower: the set of sim- 
ulated showers should by far exceed the number 
of observed showers to keep the statistical recon- 
struction errors small. 

At ultra-high energy, hybrid simulation 
schemes present a fast and efficient alternative to 
conventional Monte Carlo simulation techniques. 
In a hybrid simulation all interactions above a 
certain energy threshold are simulated with the 
Monte Carlo technique. Secondary particles that 
fall below this threshold are taken as sources of 
subshowers that are treated numerically. The 
various hybrid schemes available - for a review 
see [142) - differ mainly in the method of cal- 
culating these sub-showers. For example, the 
sub-showers could be drawn from pre-calculated 
libraries [EH] , calculated by solving cascade equa- 
tions for a Monte Carlo-generated source function 
|39j . or treated by applying numerical solutions 
of cascade equations together with analytical 
approximations or tables of shower evolution 
[143I144I14B] . 

The hybrid codes SENECA [US] and CONEX 
[39j have several interaction models implemented 
and have reached a precision and sophistication 
that makes them suited for analyzing experi- 
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mental data [TIB] , SENECA allows a full 3+1- 
dimensional simulation of EAS whereas CONEX 
is currently restricted to calculating the projec- 
tion along the shower axis. Both codes have been 
extensively compared to CORSIKA. For exam- 
ple, a comparison of the mean depth of maxi- 
mum of different models is shown in Fig. The 
agreement between the different shower simula- 
tion codes is excellent. It should also be noted 
that, only with hybrid simulation codes, showers 
can be simulated in large numbers using the time- 
consuming interaction model neXus 3.97 |147j . 

The simulation of inclined or even upward- 
going air showers has become increasingly impor- 
tant. Large- aperture experiments like Auger and 
EUSO have not only a large sensitivity to nearly 
horizontal showers but also hope to find neutrino- 
induced upward-going showers |55I148| . With 
the exception of CONEX, the currently available 
shower simulation packages are not optimized for 
such calculations. Modifications and extensions 
will be needed to allow detailed and efficient sim- 
ulation of showers at these particularly interesting 
geometries. 

4. Exotic phenomena and emulsion cham- 
ber data 

The most striking, unexpected phenomena ob- 
served in emulsion chamber experiments are Cen- 
tauro events with an exceptionally small num- 
ber of photons, events with particles or groups of 
particles being aligned along a straight line, halo 
events characterized by an unusually large area of 
darkness in the X-ray film, and deeply penetrat- 
ing cascades [14911 50) . Whether these phenom- 
ena are related to fluctuations and the measure- 
ment technique of emulsion chamber experiments 
or signs of new physics is controversially debated 
for more than 30 years. In the following we will 
briefly discuss Centauro events and comment on 
the status of searching for events with alignment 
(see |151| for a complete review). 

There are several aspects that complicate 
the interpretation of emulsion chamber measure- 
ments. First of all, many of these phenomena 
are observed only in very high energy events with 
estimated energy greater than 10 16 eV, of which 



only a small number of about 100 event has been 
collected |150| . Secondly, due to the threshold 
effect of the detectors (E 7 > 1 TeV), mainly pro- 
ton initiated events are detected |138j . As is well 
known, proton showers are characterized by very 
large shower-to-shower fluctuations (see, for ex- 
ample, |152| ). Thirdly, the number of high energy 
7-rays and hadrons cannot be measured directly - 
it is obtained by comparing the tracks at various 
depths in the detector stacks. 

Probably the most famous exotic emulsion 
chamber event is Centauro-I, an event with about 
40 high-energy hadronic jets and only one low- 
energy e.m. cluster [149] . see Fig. El Detected in 
1972 in the Chacaltaya emulsion chamber experi- 
ment |153| it represents the most extreme Cen- 
tauro event ever observed. In total there are 
about 10 Centauro events observed by the Cha- 
caltaya and Pamir experiments (see |150| . a de- 
tailed review of all events is given in |154p . No 
events of this kind were found in the experiments 
at Mount Fuji and Mount Kanbala. Also all 
searches at accelerators were negative. 

Models proposed for explaining Centauro 
events range from assuming a small fraction of 
exotic primary particles in the cosmic ray flux 
(for example, strangelets |155I156| or quark globs 
|157| ) over exotic interaction scenarios, like the 
creation of a disoriented chiral condensate 158 
or production of evaporating mini black holes by 
neutrino primaries 159 160], to conventional fea- 
tures of diffraction dissociation |161j . 

Given these ongoing attempts of explaining 
Centauro events, the results of the recent re- 
examination of the emulsion chamber plates are 
of outstanding importance |151I162| . As it turned 
out, the tracks in the two chambers that were pre- 
viously thought to belong to the Centauro-I event 
are actually due to two different events. The az- 
imuth angle of the tracks in the upper chamber 
does not match that in the lower one. As there is 
no counterpart found in the upper chamber, this 
event is still very difficult to understand |162| . 
The probability of particles produced in an inter- 
action well above the installation passing through 
a gap between the upper chambers seems to be 
very small. A point-like interaction in the upper 
layers or the wooden support frame is also mi- 
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Figure 6. Illustration of the original interpreta- 
tion of the Centauro-I event |151| . 



the KASCADE Collab. [EZj showed that aligned 
structures in hadronic shower cores at sea level 
are not related to angular correlations in hadronic 
interactions as might be expected from jet pro- 
duction 168 . Indeed the fraction of events with 
alignment is only determined by the lateral dis- 
tribution of hadrons. Measured in terms of the 
alignment parameter A4, the event distributions 
of Pamir and KASCADE data look also surpris- 
ingly similar. 

5. Gamma-ray, neutrino, and muon flux 
measurements 

Secondary particle fluxes such as hadronically 
produced gamma-rays and neutrinos provide in- 
formation on acceleration, propagation and inter- 
action of cosmic rays that are complementary to 
what can be directly deduced from the locally 
observed cosmic ray flux |169I17()I171| . In par- 
ticular, gamma-rays and neutrinos propagate on 
straight trajectories, allowing the identification of 
the source objects or environments. 



likely: the tracks would then correspond to very 
high p± particles and should point back to the in- 
teraction vertex. No such geometric convergence 
of the tracks is found. This means that many par- 
ticles of almost parallel trajectories have hit the 
emulsion chamber without interacting in the up- 
per lead stack, again a very exotic scenario that 
lacks explanation. 

The situation is similarly controversial regard- 
ing the experimental information on events with 
co-planar particle emission. 

The Pamir Collab. find an excess of events 
with alignment of the substructures for E > 8 — 
10 x 10 15 eV |163I164| . The highest energy event 
measured with an emulsion experiment during a 
series of Concorde flights (average atm. depth 
~ 100g/cm 2 ) shows also impressive alignment 
165 . At lower energy, no excess of coplanar par- 
ticle production is found. For example, measure- 
ments in the energy range below 10 14 eV by the 
RUNJOB Collab. |166j and also a direct search 
with the CERN NA22 experiment at 2.5 x 10 11 eV 
|164j provided distributions that agreed with the 
expectations. Furthermore, a recent study of 



5.1. Gamma-rays 

With the begin of routine operation of 
the second generation imaging atmospheric 
Cherenkov telescopes (IACT) CANGA- 
ROO [IZ21, HESS |17SI174| . MAGIC |l75ll76j . 
and VERITAS |177| . many new TeV gamma-ray 
sources are discovered and their spectra mea- 
sured. It is impossible to summarize the progress 
in this extremely active and diverse field of re- 
search and any comments made here will be soon 
outdated. 

At the time of writing this article all big four 
IACT installations arc completed and take data. 
Whereas CANGAROO and HESS have already 
several telescopes online, MAGIC and VERITAS 
work with single, but bigger telescopes. Both the 
MAGIC and VERITAS Collaborations are in the 
process of adding another telescope for stereo- 
scopic observation with greatly improved back- 
ground rejection. The HESS telescope system is 
characterized by an unprecedented angular reso- 
lution of 0.06°. The MAGIC telescope has a light- 
weight design for very fast slew to observe tran- 
sient sources. It's low-energy threshold is planned 
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to be about 20 GeV as compared to ~ 70 GeV for 
HESS and VERITAS 175 . The four telescopes 
together give almost uniform full-sky coverage. 

Highlights of the early HESS data taking are 
certainly the measurement of the gamma-ray flux 
from the Galactic Center |178j and the first obser- 
vation of a SNR as spatially resolved TeV gamma- 
ray source (RXJ 1713-3946, a possible site of cos- 
mic ray acceleration) |179| . Both sources were 
previously detected with the CANGAROO tele- 
scopes 180 l HT] but with much more limited res- 
olution. The potential of the HESS telescopes is 
also underlined by the serendipituous discovery of 
an unknown TeV gamma-ray source, now called 
TeV J1303-63, in the field of view of the binary 
pulsar system PSR B1259-630 |T75| . 

In contrast to imaging Cherenkov telescopes, 
air shower arrays can be used to continuously 
monitor the gamma-ray sky with very high duty 
cycle and wide field of view. The two cur- 
rently operated installations of this type are Ti- 
bet AS 7 |22ll82| and Milagro |183ll84| being lo- 
cated at an altitude of 4300 m and 2350 m, re- 
spectively. Whereas the Milagro Collab. employ 
an active hadron/gamma-ray separation via two 
layers of PMTs in an 8 m deep water pond, the 
Tibet experiment searches for arrival direction 
anisotropies due to gamma-rays on top of the 
isotropic cosmic ray background with a dense 
scintillator array. 

Both experiments have detected the Crab SNR 
and the active galaxy Mrk 421 at the 5a level 
[1851186) . New results from Milagro are the de- 
tection of TeV gamma-rays from the entire inner 
Galactic plane region and the observation of two 
extended sources, one coincident with EGRET 
source 3EG J0520+2556 and another one in the 
Cygnus region of the Galactic plane |187I188| . 
It is intriguing that the latter source coincides 
within the Milagro resolution of about 2° with 
the HEGRA source TeV J2032+4131 HSU and 
the region from where AGASA reported an excess 
of ~ 10 18 eV cosmic rays 190 . At the moment 
an interpretation of these observations in terms 
of a very high energy cosmic ray source (region) 
is too speculative - more data will be needed. 



5.2. High-energy neutrinos 

The interpretation of gamma-ray fluxes from 
potential cosmic ray sources suffers from am- 
biguities due to the superposition of different 
gamma-ray production processes: ir° decay, in- 
verse Compton scattering, synchrotron radiation, 
and bremsstrahlung. These uncertainties are ex- 
pected to be much smaller for neutrinos as they 
are mainly produced in hadronic interactions via 
the decay of pions and kaons. Furthermore neu- 
trinos can travel over large distances virtually 
unattenuated and are, therefore, ideal messenger 
particles, allowing a multitude of astrophysical in- 
vestigations |191| . On the other hand, their small 
interaction cross section requires very large effec- 
tive detector volumes. 

There are two neutrino telescopes taking data 
at the moment, AMANDA-II |62I192| and Baikal 
NT-200 |193ll94j . Although limited by detec- 
tor size, the sensitivity of both detectors will ap- 
proach the cascade bound in the next years 3 , 
i.e. touch the region where one can hope for 
a discovery. The Waxman-Bahcall bound 196 
(see also discussion in |197| V often considered 
as a reference flux that is guaranteed if pro- 
tons are the ultra-high energy cosmic rays, can- 
not be reached with these installations. About 
3300 (370) neutrino candidates are found in the 
AMANDA (Baikal) data taken until end of 2003. 
The number of neutrinos and their distribution 
are compatible with the expectation from atmo- 
spheric neutrino production - neither a diffuse 
flux of extra-terrestial neutrinos nor astrophysi- 
cal point sources have yet been discovered. 

The construction of the successor to AMANDA 
and much bigger neutrino telescope IceCube is on 
track |198| . IceCube will have a sensitivity that 
reaches well into the region below the Waxman- 
Bahcall bound, promising discoveries and many 
astrophysical and particle physics applications 
|199| . The first IceCube string was successfully 
deployed in February 2004 P%| . 

The Mediterranean neutrino telescope col- 
laborations (ANTARES [230], NESTOR jM], 



3 The cascade bound, also called gamma-ray bound, is 
based on the assumption that all observed extragalac- 
tic gamma-rays were produced together with neutrinos in 
hadronic cascades 117111951 . 
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NEMO 202 1) [T9"3) have performed prototype in- 
stallations and test runs. In 2003 the ANTARES 
Collab. operated a prototype sector line with 
PMTs and a mini instrumentation line at the se- 
lected ANTARES site near Toulon (2500m water 
depth) . Valuable information on the performance 
of cables and connectors under the harsh deep-sea 
conditions were gathered. It is planned to build 
the complete ANTARES detector of 12 strings 
and in total 900 PMTs in 2005 - 2007. NEMO has 
selected a site close to Sicily with nearly perfect 
conditions at 3500m water depth which is contin- 
uously monitored. The project is in the advanced 
R&D stage with the plan to build a prototype in 
2005. The NESTOR site provides a large plateau 
at the sea floor at about 4000 m water depth. 
In 2003 one fully equipped prototype "star" of 
32 m diameter of a NESTOR tower was success- 
fully operated, allowing the measurement of the 
atmospheric muon flux. It is planned to install 7 
complete towers by the end of 2006, providing a 
detector of about 0.15 km 3 . 

It is clear that a km 3 -sized neutrino detector 
is needed in the northern hemisphere to comple- 
ment the field of view of IceCube. Therefore the 
Mediterranean neutrino telescope collaborations 
have recently joined their efforts to construct such 
a detector by initiating the design study KM3NeT 

To measure neutrino fluxes at even higher en- 
ergy, radio emission of neutrino-induced showers 
in dense materials can be employed 204 . Several 
experiments have recently performed measure- 
ments and derived first limits on the neutrino flux 
at ultra-high energy (FORTE [gfiS] . GLUE |2JE|, 
ANITA For example, ANITA is designed 

to search for neutrinos with E u > 3 x 10 18 eV 
by monitoring radio signals from the antarctic 
ice cap using a balloon-born system of antennas 
204 207 . A preparatory test flight with a pro- 
totype instrument (ANITA-lite) was performed 
during the 03/04 austral season [204) . Already 
on the basis of the prototype flight, giving about 
7 days of data, a competitive limit on the ultra- 
high energy neutrino flux could be derived 204 . 



5.3. High-energy muons 

Atmospheric muons, being a major background 
for neutrino telescopes, carry valuable informa- 
tion as messengers of hadronic interactions in the 
atmosphere. Muons are also directly linked to at- 
mospheric neutrino production and can be used 
to test predictions on neutrino fluxes as needed 
for oscillation parameter analyses [2088209) . Par- 
ticularly interesting is the comparison of muon 
fluxes to simulations performed with the same 
codes as used for air shower analyses [210121112091 
1212) . Of course, muon flux predictions depend on 
both the used hadronic interaction model and the 
assumed primary cosmic ray flux. 
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Figure 7. Vertical atmospheric muon flux as mea- 
sured by L3+C. The upper panel shows the flux 
of all muons and the lower panel the charge ratio. 
The data are compared with different theoretical 
predictions [2101213) . 



Very precise high-energy muon measurements 
can be carried out by particle physics detectors 
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at colliders. For example, Fig. [7| shows the in- 
clusive muon flux and charge ratio measured by 
the L3 Collab. |213| . The experimental results for 
vertical muons are compared to different theoret- 
ical predictions. In this case, in contrast to EAS 
simulations, QGSJET01 predicts a smaller muon 
flux than SIBYLL 2.1 (cmp. Fig. EJ). Hadronic 
interaction models tuned for muon and neutrino 
flux calculations give a better description of the 
data |212| . As known from simulations at lower 
energy |2()HI214| . the muon charge ratio is found 
to be very sensitive to the production of fast 
secondary hadrons. None of the models imple- 
mented in CORSIKA gives a good overall descrip- 
tion of the data. The deficit of muons found rela- 
tive to CORSIKA simulations is also seen by two 
other LEP experiments. The DELPHI and Alcph 
groups find that the number of high energy muon 
bundles cannot be described by the hadronic in- 
teraction models available in CORSIKA, assum- 
ing even a completely iron-dominated primary 
composition |215I216| . 

At higher energy, the AMANDA and NESTOR 
collaborations have also measured the atmo- 
spheric muon flux in the TeV energy region 
62 153]. These measurements are integral flux 
determinations because of the very limited energy 
resolution of neutrino telescopes. The AMANDA 
Collab. have also compared their measurement to 
simulations with CORSIKA and find good agree- 
ment within the experimental uncertainties |62| . 

6. Conclusions and outlook 

Flux and composition of ultra-high energy cos- 
mic rays are still very uncertain because of the 
low statistics of showers observed so far and the 
model-dependence of the shower data interpre- 
tation. Significant progress in this field is ex- 
pected by new large-aperture installations - the 
Pierre Auger Observatory and the Telescope Ar- 
ray [SSj. To study the flux at even higher en- 
ergies ~ 10 21 eV with sufficient statistics, new 
techniques will be required. Observing the atmo- 
sphere from outer space is one possible solution to 
increase the aperture further (for example, EUSO 
|148| 1. Another possibility could be the use of 
radio antenna arrays similar to that of particle 



detectors PT7] . 

The situation is similar in the knee energy re- 
gion. Although the all-particle flux is known 
rather well there are large uncertainties in the 
composition. Nevertheless a clear trend from a 
mixed composition at low energy to a predom- 
inantly heavy one above the knee is seen in all 
recent measurements. The experimental errors 
are completely dominated by the systematic un- 
certainties due to our limited understanding of 
hadronic interactions at high energy and in for- 
ward direction. 

The dependence on air shower simulations sim- 
ulations can be reduced by combining different 
detection techniques to measure qualitatively dif- 
ferent shower observables at the same time. Ex- 
periments of this type are the Pierre Auger Ob- 
servatory, TA, and IceCube/IceTop. 

Whereas there are several detectors measuring 
the cosmic ray flux in the knee region and at 
ultra-high energy, there is a lack of data in the 
energy range 10 17 — 10 19 eV. It is clear that the 
latter range is of great interest as it expected to 
cover the transition from Galactic to extragalactic 
cosmic rays. KASCADE-Grande and IceTop will 
measure showers only up to 10 175 eV with good 
statistics. Therefore it is worthwhile to upgrade 
large-aperture instruments such as Auger or TA 
to extend their energy range down to 10 17 eV. 

The field of emulsion chamber measurements 
is still full of mysteries. After more than 30 
years the interpretation of one of the most fa- 
mous emulsion chamber events, Centauro-I, has 
changed completely, now being even more diffi- 
cult to explain. Not a single non-emulsion ex- 
periment could confirm any of the claimed exotic 
event features. Substantial progress in this field 
can only be expected by new measurements com- 
bining large- aperture emulsion stacks and modern 
particle detectors. 

Many of the questions related to cosmic rays 
and astroparticle physics can only be solved by 
measuring and understanding secondary particle 
fluxes. There has already been enormous progress 
in the field of gamma-ray and neutrino measure- 
ments and much more can be expected for the 
next years. The second generation imaging air 
Cherenkov telescopes will provide high resolution 
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images of TeV gamma-ray sources and water/ice 
detectors of km 3 size will probe the neutrino flux 
in the same energy range. Neutrino fluxes at 
ultra-high energy will be searched for by large- 
aperture air shower installations and dedicated 
radio signal experiments. The field of high-energy 
neutrino astronomy is still in an infantile stage 
but with a bright future ahead. 

For all these measurements and related data 
analyses, detailed simulations are very impor- 
tant. Shower simulation tools in general and 
hadronic interaction models in particular should 
be improved continuously and tested by compar- 
ing them with a large variety of data. During the 
last decade great progress was achieved by the 
introduction of multi-purpose code packages such 
as CORSIKA and AIRES that are professionally 
maintained. However, it should not be overlooked 
that the quality of air shower and inclusive flux 
simulations depends crucially on particle produc- 
tion data measured in fixed-target and collider 
experiments. At the moment the lack of suit- 
able accelerator data is the dominant source of 
systematic uncertainties in cosmic ray measure- 
ments. As we don't have a calculable theory 
of hadronic multiparticle production, there is no 
change of this dependence on accelerator data to 
be expected in near future. 
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